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Abstract 

The electric field generated by a time varying point charge in a three-dimensional, unbounded, spatially 
homogeneous plasma with a uniform background magnetic field and a uniform (static) flow velocity is stud- 
ied in the electrostatic approximation which is often valid in the near field. For plasmas characterized by 
Maxwell distribution functions with isotropic temperatures, the linearized Vlasov-Poisson equations may 
be formulated in terms of an equivalent integral equation in the time domain. The kernel of the integral 
equation has a relatively simple mathematical form consisting of elementary functions such as exponential 
and trigonometric functions (sines and cosines), and contains no infinite sums of Bessel functions. Con- 
sequently, the integral equation is amenable to numerical solutions and may be useful for the study of the 
impulse response of magnetized plasmas and, more generally, the response to arbitrary waveforms. 
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I. INTRODUCTION 



The study of plasma waves generated by an oscillating point charge or point dipole — fixed in 
space but oscillating in time — in a spatially homogeneous hot magnetized plasma is an important 
model problem that has been investigated both theoretically and experimentally since the 1960's 
1 1-39]. Following the pioneering work of Landau in 1946 [40], the fully kinetic self-consistent 
field approach described by the coupled Vlasov and Maxwell equations has been used to solve 
and study such systems, mainly through Fourier and Laplace transform analysis. Experience with 
these techniques shows that even in the simplest spatially homogeneous systems, the presence of 
ambient magnetic fields usually complicates the mathematical analysis significantly. 

Here we show that in the electrostatic approximation and for particle distribution functions 
with isotropic temperatures, an alternative approach based on solving an integral equation in the 
time domain provides a relatively simple mathematical formulation that contains no Bessel func- 
tion series. Similar formulations may be possible for gyrotropic distributions with anisotropic 
temperatures, but these are not considered here. To the author's knowledge, the integral equation 
presented here for the field of a point source in a magnetized plasma has not appeared before in 
the literature. 

Solutions of the linearized Vlasov-Maxwell or Vlasov-Poisson equations for specified charge 
and current distributions are of interest in many applications where the steady state solution for 
time harmonic forcing and the impulse response corresponding to a delta function forcing (in time) 
are usually of special interest. In the electrostatic approximation and in the presence of a constant 
background magnetic field, the steady state solution for a point source can usually be expressed 



in terms of an infinite series of modified Bessel functions [119 



30 



34 



4ll44-3ll. Here, it is shown 



that the solutions for a point source with any prescribed time dependence can be derived from 
an integral equation in which the kernel is free of such Bessel function series. Since the integral 
equation formulation is amenable to numerical calculations, it can be used to study the impulse 
response of magnetized plasmas as well as the response to other arbitrary waveforms through 
direct calculation. In that respect, it has more flexibility and some other advantages compared to 
the Laplace transform approach. 
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II. STATEMENT OF THE PROBLEM 



Consider a spatially unbounded and homogeneous plasma with a background magnetic field 
B = B e z , B > 0, and an externally imposed point charge with charge density q(t)5(x), where 
q(t) is a given function of time, 5(x) = 5(x)8(y)5(z) is the three dimensional delta function, 
8(x) is the (one dimensional) Dirac delta function, and (x, y, z) are the usual orthogonal cartesian 
coordinates in three dimensional space with corresponding unit vectors e x ,e y ,e z . The function 
q(t) is assumed to vanish for t < and is continuous for t > 0. Hence, by causality, the plasma 
response will vanish at t — or, in other words, the initial conditions on the perturbed fields and 
the perturbed distribution functions all vanish at t = 0. For t > 0, the point charge excites plasma 
waves which propagate away from the origin and the problem is to compute the fields produced 
by this point source. 

At equilibrium, the plasma is charge neutral and current free which is expressed by the relations 

^2n s q s = and ^n s q s V s = (1) 

s s 

where n s , q s and V s are the equilibrium number density, charge, and bulk flow velocity of particle 
species s, respectively. The equilibrium distribution functions are all assumed to be non-relativistic 
convected Maxwell distributions with isotropic temperatures of the form 

f° s ( v } = 7 — ^I7^ ex P 

where v s = (2k B T s /m s ) 1 ^ 2 is the thermal speed, k B is Boltzmann's constant, m s is the particle 
mass, and T s is the kinetic temperature. Note, however, that the analysis below can readily be 
generalized to other distribution functions with isotropic temperatures. The Maxwell distribution 
© has the convenient property 

V v f 0s = -^(v-V s )f 0s (v) (3) 

that shall be used below. 

To linearize the Vlasov-Poisson equations, the first step is to specify the equilibrium state. In 
the presence of plasma flow in an arbitrary direction, the equilibrium solution of the collisionless 
Vlasov equation must satisfy (E + v x B ) ■ (v — V s ) = which implies the existence of a static 
electric field E = —V s xB . For the solution to be self consistent, the component of the velocity 
transverse to B must be the same for each species, that is, V s x B must be the same for each 
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(2) 



s. The equilibrium state consisting of E , B , and fo(v) can now be used to linearize the Vlasov 
equation. To simplify the presentation, the ion response shall be neglected from now on and only 
the electron contribution shall be considered; because the susceptibilities of the different particle 
species are additive, the ion response can easily be taken into account as shown at the end of the 
derivation. 



The perturbed electric field E\ is assumed to be electrostatic meaning that E\ 



yd> , w 



here 



12b, at 



(j)(x, t) is the electrostatic potential, an approximation that is valid in the near field [|12l 
least over some range of physical parameters. The precise range of validity of the electrostatic 
approximation for the problem under consideration is unknown and shall not be studied here. 
However, for time harmonic forcing, the steady state solution can be computed numerically both 
with and without the electrostatic approximation, thus providing a means of comparison. 

In terms of the velocity variable u = v — V , where the subscript on V = V e is omitted for 
convenience, the linearized Vlasov equation takes the form 

^ + (« + V) ■ VA - — (u x B ) ■ V u /i = -— V0 • V u /o, (4) 
at m e m e 

where fi(x, u, t) is the perturbed distribution function defined such that the complete distribution 
function is n e (f + fi) with fo(u) normalized to unity and e > the electronic charge. Express 
the velocity vector in cyclindrical coordinates as 

u = u± cos((p)e x + u± sm(<p)e y + u\\e z . (5) 

After Fourier transformation with respect to the spatial variables, the Vlasov equation becomes 



— + ik ■ [u + V)fi - ^ e ^— = -5 — ik ■ uf 



... » • ..i i --. . ■- ~Jo{u), (6) 

at otp vpn e 

where fi(k, u, t) and (p(k, t) are the respective Fourier transforms of fi(x, u, t) and (f)(x, t), and 
VL e = —eBo/m e is the signed electron cyclotron frequency (fl e < 0). The problem is to solve the 
Vlasov equation (J6]) for t > together with Poisson's equation 

k ^(k,t) = -^- [ f 1 (k,u,t)du + ^- (7) 

^0 J e 

subject to the initial conditions fi(k,u,0) = and 0(fe, 0) = 0, where q(t) is a continuous 
function such that q(0) = 0. Otherwise, the forcing function q(t) is arbitrary. 
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III. FORMULATION AS AN INTEGRAL EQUATION 



The goal is to reduce the initial value problem consisting of the Vlasov-Poisson equations © 
and © to an equivalent integral equation. The result is as follows. The potential (f>(k, t) which 
solves the Vlasov-Poisson system is also the solution of the Volterra integral equation of the second 
kind 



k,t)+u 2 T{k,t-T)<j){k,T)dT 



9C0 



MO, 



(8) 



with kernel 



r(fc,r) 



sm(f2 e r) k 2 



x exp < — 



k\v 2 



e sin 2 



knVfT 



ik ■ Vt 



(9) 



Note that the kernel is a convolution kernel as is required for a linear time invariant system. Re- 
markably, the kernel © does not contain the infinite sums of Bessel functions that occur in steady 
state solutions with time harmonic forcing. The terms containing fi e represent the magnetic field 
effects. In the limit as Vl e — > 0, the function © reduces to the correct kernel for a point charge in 
an unmagnetized plasma 

2 



r(fc, r) = r exp 



kv e T 



ik ■ Vt 



(10) 



where k = (k\ + k 2 ) 1 / 2 . Because the kernel © is smooth (continuously differentiable) and con- 
sists of compositions of algebraic, exponential and trigonometric functions (sines and cosines), 
it is convenient for purposes of numerical calculation. When only the electron response is con- 
sidered there are two inherent timescales which must be resolved, the electron cyclotron period 
2n/co ce and the inverse plasma frequency 2n/oj pe . 

Integral equation formulations of the Vlasov-Poisson system in unmagnetized plasmas are well 
known Il44l - l52ll . However, to the authors' knowledge, the integral equation formulation for a 
magnetized plasma given by ([8]) and © is new. Of course, the literature is vast, spanning more 
than 50 years, and it is possible that this result may have been derived previously in work unknown 
to the present author. The remainder of this paper is devoted to a derivation of the integral equation 
© and the kernel ©. 
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IV. SOLUTION OF THE VLASOV EQUATION 



Similar to ©, let 

k = k± cos(9)e x + kj_ sm(6)e y + kne z . (11) 

Then 

k ■ u = kj_u± cos(f — 9) + k\\U\\. (12) 

The substitution 

/i = ^exp| + sm(y> - 0) -i{k\\ u \\ + k ' ( 13 ) 

brings the Vlasov equation © into the form 

dip dip 2c f k j_Hj _\_ 1 

Consider the transformation from independent variables t and ip to a new set of independent vari- 
ables 77 and £ defined by 

77 = |0 + fi e t) 



It follows that 



(15) 



dip dip di) 

^~ n ^ = - n ^ (16) 



and, therefore, the solution of (jT4b is obtained by integrating both sides of that equation with 
respect to £ while holding 77 fixed. Because t = is equivalent to £ = 77, the integration is 
performed from the lower limit £' = 77 to £' = £, where £' is the variable of integration. To express 
the right-hand side of (fT4l) in terms of 77 and £, substitute ip = 77 + £ and i = (77 — £)/fi e so, for 
example, 

k ■ u = k±u± cos(?7 + £ — 9) + &||it||. (17) 
The integration of ([14)) implies 

2e ftdt'~ (l _ 77 

-lU / , f/ M , I 1 / ■ Sm (^ + £' ~ A) + fcp|| (g ~j) ) , 1fi . 

x 7 cos(77 + £ — 0) + «||W||J exp < — ? — — >. (18) 
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By definition, ipij], rj) is equal to ip(t = 0) which is equal to zero by virtue of the vanishing initial 
condition on f\. The change of variable r = (77 — £')/f2 e , where rj is a constant, yields 

2e ~ 

exp I - z sin(2?7 - - fi e r) + i(Jfe|| U || + k ■ V)r J. (19) 

Because 77 is held constant inside the integrand, if one makes the substitution r] = (ip + Cl e t)/2 
inside the integrand with both ip and t held constant during the integration, then the result of the 
integration is the same as if the substitution 77 = (y? + Vt e t)/2 were made after the integration. 
Therefore, the preceding equation is equivalent to 

2e f* ~ 

ip(k, u, t) = — / dr </>(k, T)f (u)i{ k_i_u± coshp — + QJt — r)l + k\\U\\ } 

vjm e J 

exp I - i sin[(p - + Q e (t - r)] + i{k\\U\\ + k ■ V)r J. (20) 

It can be verified by direct calculation that this is the solution of the partial differential equation 
(PT4l) . Thus, from (fT3"T ), the solution of the Vlasov equation © is 

~ 2e f f ~ 

fi(k, u, t) = — — / dr <p(k, t) f (u)i{k±u± cos ho — + VtJt — r)] + k\\U\\\ 
vim e J 



exp 



L e JO 

k±u± r . 



or, equivalently, 



I sinks -0 + Q e (t - t)] - sin(^ - 0)} - i{k\\U\\ + k ■ V)(t - r) 



(21) 



~ 2e /' ~ 

fi(k,u, t) = — — / dr </>(k, t — r)f (u)i \k±u± cos(y? — + Q e r) + k\\u\\\ 
vjm e Jo 

\ Jv _\_U _\_ I 

exp < — i — — — [sin(y? — 8 + Vt e r) — sin(<y9 — 0)] — i(k\\u\\ + k ■ V)r > (22) 
This result is now inserted into Poisson's equation © and the integration over velocity space is 



performed. 



V. VELOCITY SPACE INTEGRALS 



The integration of / x with respect to u is performed in cylindrical coordinates. To facilitate the 
calculation, it is expedient to write ([22)) in the form 



f 1 (k,u,t) = --^- ! dr<P(k,t-r)f (u)e- ik - VT 
v 2 e m e J 

< — i — — — [sin(y? — + fi e r) — sin(y? — 0)] — zfe[|U||T > (23) 



x ^-exp 
or 
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The integal with respect to <p is readily computed using the Bessel function series 



J2 e ^ J - 



e lzsm <>= > e--J n iz) 



(24) 



where J n (z) is the Bessel function of the first kind of order n and argument z, and z is an arbitrary 
complex number y53|]. It follows from this series expansion that 



exp 



i srnO - 9) 



and 



cxp 



i sin(y? — 9 + fl e r) 



oo , 

c mfy-g+n e r) ^ fc -L^-L ^ 



(25) 



(26) 



The orthogonality relations for the complex exponentials then yield 



o 



2 \ S 

dip exp < — % — — — 

I i 



[ sm((p — 9 + fi e r) — sin(y9 — 0)] j- 



/l — — oo 

The series on the right-hand side is a special case of the Fourier series 



(27) 



2 = J n (2zsin 



E ^ 

71= — OO 

with Fourier coefficients given by I54I page 48, equation 13] 



(28) 



- / Jo^sin-le 



-in<j) 



[Jn 



(29) 



Hence, 

f27T 



dy> exp |-i — — — [ sin(<£> — 6 1 + fi e r) — sin(<£> — #)] j 



2^J [-—sm — 



(30) 



and, from (1231 . 

»2tt 



d ¥ ,/ 1 (fe,«,t) = _-_ / dr0(fe,t-r)/ o («)e-* w 
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It remains to perform the integration with respect to u± and u\\. The integration with respect to 
u± is accomplished by means of Weber's first exponential integral J55, page 393] 

— 2 \ ± e<rtjo(au ± )du ± = ^exp( - ^-X. (32) 
nvlJo 2tt [ 4 J 

The application of (l32l) yields 

-2tt r. /-t 



/ u±du ± dtp fi(k,u,t) = — - — / dr(j)(k } t-T)f (u\\)e 
Jo Jo v e m e J 



—ik-Vr 



And the integration with respect to u\\ now gives 



/ u±du± / ditii / d(p fi(k,u,t) = — - — / dr (f)(k,t — r)e~ lk l 7 
Jo ' J -00 Jo v 2 e m e J 



d ( k\vl o {Q p t\ fk\\v P r\ 2 
^ 6XP 



x _exp<!-^sin 2 ^ - P^-l \. (34) 



Differentiation with respect to r yields the final result 

r _ ^2 e rt 

f 1 (k,u,t)du= — / drT(k,r)(j)(k,t-T), (35) 
J rn e J 

where T(k,r) is the kernel ©. Substitution of this result into Poisson's equation © yields the 
integral equation dH). This completes the derivation of dE) and ©. 



VI. INTEGRAL EQUATION FOR MULTI-SPECIES PLASMAS 

Different types of ions and other charged particle species can be taken into account in a straight- 
forward manner. The Vlasov equation for particle species s takes the form 

df\ ~ df\ 2q (f) 

— - +ik- (u + V s )fu - ^ = — ^—ik ■ uf 0s (u), (36) 
at 0(p vjm s 

where f ls (k,u,t) is the Fourier transform of f ls (x,u,t) and f2 s = q s B /m s is the signed cy- 
clotron frequency of species s. The effects of the different species are additive in Poisson's equa- 
tion which becomes 

*^E- / fis(k,u,t)du + ^l. (37) 
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Subject to the initial conditions f ls (k,u,0) = and </>(&, 0) = 0, the Vlasov equation (1361) is 
solved as described in Section [IV] and the velocity space integrals are computed as in Section |V] 
with the result 



where 



r(fe, t — r)(j)(k, t) dr 



(38) 



sm(Q s r) k\ 
SI,t k 2 



x exp < — 



kjv? 



sm 



k\\V<,T 



ik-V s r\ (39) 



and u 2 = n s q 2 /eom s . Even though the integral equation formulation of the Vlasov-Poisson 
system is exact, the calculation of numerical solutions in multi-species plasmas is difficult in 
practice because of the disparity between ion and electron timescales. 



VII. CONCLUSIONS 

It has been shown that the electric field generated by a time varying point charge q(t)S(x) in 
a magnetized plasma can be obtained as the solution of an integral equation and, in the case of a 
convecting Maxwellian plasma with an isotropic pressure distribution, that the kernel may be ex- 
pressed in terms of elementary functions. Numerical solutions of this integral equation can be used 
to study the response of the plasma to various types of forcing q(t) and thus it provides an alterna- 
tive approach to exact analytic theories and their approximations which are often cumbersome for 
the analysis of magnetized plasmas. For example, the impulse response of a magnetized plasma 
can be studied by using pulsed sources such as q(t) = Atexp(—t/to) or q{t) = Atexp(—t 2 /2tl). 
These functions both rise to a maximum at t — to and then decline to zero; the risetime to is ad- 
justable, of course, for the application of interest. The integral equation can also be used to study 
wave excitation and propagation for sources with arbitrary waveforms such that q(t) is continuous 
and satisfies q(0) = 0. The integral equation, therefore, should be of practical value. 

Concerning the avoidance of infinite series of Bessel functions in studies of homogeneous mag- 



netized plasmas, the recent work by Qin et al. pa] should be mentioned. Textbook derivations of 



the electromagnetic susceptibility of homogeneous, magnetized, gyrotropic plasmas usually em- 
ploy expansions in infinite series of Bessel functions to perform the integration over the particle 
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gyro-phase [see, for example, chapter 10 in Ref. 57J. Qin et al. Il56h have shown how this integra 



tion may be carried out without the use of Bessel function expansions. While this is an important 
mathematical development, the resulting form of the susceptibility tensor derived by Qin et al. 
still contains integrations over the variables p± and p» which have not been carried out. In the case 
of Maxwell or bi-Maxwell velocity distribution functions, the remaining integrations over p± and 
p» inevitably give rise to infinite series of Bessel functions and lead to the well known expressions 



given, for example, by Stix [15711 . Consequently, for Maxwellian plasmas, the expressions for the 
susceptibility tensor derived by Qin et al. have postponed but not eliminated the practical need for 
such series expansions in numerical calculations of the hot plasma dielectric tensor. 
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